WHERE DO COOLING FLOWS COOL? 
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ABSTRACT 

Typically 5 percent of the total baryonic mass in luminous elliptical galaxies is in the form 
of cooled interstellar gas. Although the mass contributed by cooled gas is small relative to the 
mass of the old stellar system in these galaxies, it is almost certainly concentrated within the 
optical effective radius where it can influence the local dynamical mass. However, the mass of 
cooled gas cannot be confined to very small galactic radii (r < O.Olre) since its mass would 
greatly exceed that of known central mass concentrations in giant ellipticals, normally attributed 
to massive black holes. We explore the proposition that a population of very low mass, optically 
dark stars is created from the cooled gas. For a wide variety of assumed radial distributions for the 
interstellar cooling, we find that the mass of cooled gas contributes significantly (^^30 percent) to 
stellar dynamical mass to light ratios which, as a result, are expected to vary with galactic radius. 
However, if the stars formed from cooled interstellar gas are optically luminous, their pertubation 
on the the mass to light ratio of the old stellar population may be reduced. Cooling mass dropout 
also perturbs the local apparent X-ray surface brightness distribution, often in a positive sense 
for centrally concentrated cooling. In general the computed X-ray surface brightness exceeds 
observed values within r^, suggesting the presence of additional support by magnetic stresses or 
non-thermal pressure. The mass of cooled gas inside r^ is sensitive to rate that old stars lose mass 
Af*, but this rate is nearly independent of the initial mass function of the old stellar population. 



Subject headings: galaxies: elliptical and lenticular 
galaxies: interstellar medium - x-rays: galaxies 

1. INTRODUCTION 

Perhaps the most perplexing and long-standing 
problem associated with galactic and cluster cool- 
ing flows is the uncertain physical nature and spa- 
tial distribution of the gas that cools. The ap- 
parent absence of large masses of cooled gas in 
elliptical galaxies has led some to argue that lit- 
tle or no cooling actually occurs and to postulate 
some source of heating that offsets the radiative 
losses in X-ray emission. But the energy required 
to balance radiative losses is prohibitively large 
and appropriate heating sources may not be uni- 
versally available. If the expected radiative cool- 
ing actually occurs, two questions arise: (1) What 
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is the nature of the objects that condense from 
the hot gas? and (2) Where is most of the cooled 
mass located in the galaxy? Regarding the first 
question, a variety of physical arguments support 
the hypothesis, or even the inevitability, of low 
mass star formation (Fabian, Nulsen & Canizares 
1982; Thomas 1986; Cowie & Binney 1988; Ved- 
der, Trester, & Canizares 1988; Sarazin & Ashe 
1989; Ferland, Fabian & Johnstone 1994; Mathews 
& Brighenti 1999). Here we shall address the sec- 
ond question in the context of cooling flows in el- 
liptical galaxies where the known stellar mass and 
light profiles strongly constrain the spatial distri- 
bution of cooled gas. 

We adopt the generally accepted hypothesis 
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that only stars of very low mass (e.g. ^0.1 Mq) 
form in cooling flows (e.g. Ferland, Fabian & 
Johnstone 1994), so that the mass to light ratio 
of the young stellar population formed from the 
cooling gas is essentially infinite. In view of the 
difficulties we encounter with this hypothesis, de- 
scribed below, it seems more likely that the stel- 
lar population formed from cooled gas extends to 
somewhat more massive stars that are optically 
luminous. 

Our gas dynamical models for the evolution of 
hot interstellar gas in giant ellipticals indicate that 
the origin of the gas varies with galactic radius. 
Most of the gas in the inner, optically luminous 
regions originates from the ejected envelopes of 
evolving stars; gas in the outer halo is supplied by 
cosmological secondary infall or tidal acquisitions 
from neighboring galaxies (Mathews & Brighcnti 
1998b). Circumgalactic gas around massive ellip- 
ticals is enriched by Type II supernovae that ac- 
companied early star formation. The variability 
of circumgalactic gas among luminous ellipticals 
is responsible for some of the enormous dispersion 
in X-ray luminosity among ellipticals of simi- 
lar optical luminosity Lb (Mathews & Brighenti 
1998a). 

Since the hot interstellar gas in a bright ellip- 
tical emits observable X-rays, it is clearly losing 

energy. However, as the gas loses energy it is com- 
pressed toward the galactic center by gravitational 
forces and Pdv work maintains the high temper- 
atures observed, T ^ 10^ K, producing a galactic 
cooling flow. The positive interstellar temperature 
gradients typically observed within a few effective 
radii are often cited as evidence of radiative cool- 
ing in a cooling flow, but this cooling is due instead 
to the mixing of somewhat cooler, locally virial- 
ized gas ejected from stars with hotter gas arriving 
from larger galactic radii (Mathews & Brighenti 
1998b; Brighenti & Mathews 1998, 1999a). If 
large entropy fluctuations are present in the hot 
gas, catastrophic cooling can occur at any radius 
in the flow. Regions of low entropy (low tempera- 
ture, high density) radiate more and cool sooner. 
The amplitude distribution of entropy fluctuations 
in the interstellar gas determines the radius where 
cooling mass dropout occurs in the cooling flow. 
For example, if the entropy in some region in the 
flow is only slightly less than in the ambient flow, 
the differential radiative cooling will be small and 



the region will cool out of the flow at small galac- 
tic radii; conversely, localized regions with entropy 
much lower than the ambient flow cool rapidly and 
deposit their mass at large radii. Some possible 
sources of interstellar entropy variations are stel- 
lar winds, explosions of Type la supernovae, non- 
uniform SNII heating at early times, and mergers 
with small, gas-rich galaxies. 

The total rate that mass cools and drops out of 
the flow is closely related to the X-ray luminosity 
Lx- The X-ray luminosity can be approximately 
expressed as the product of the total cooling rate 
M and the enthalpy per gram in the hot gas, or 

M={j^) i.,bo; «^ 2.5Meyr-\ 

Here wc have used data from the giant Virgo el- 
liptical NGC 4472: T w 1.3 x 10^ K; L^(0.5 - 
4.5 keV) = 4.5 x 10*^ L^^boi « 1.6L^(0.5 - 
4.5 keV). If Lx and T are reasonably constant 
over the Hubble time, a mass Meg ~ 3 x 10^° Mq 
of cold gas is expected to condense from the hot 
ISM somewhere within NGC 4472. Although this 
mass is very large, it is only about 4 percent of the 
total stellar mass in NGC 4472 today. The mass 
that cools can therefore be ignored if it is widely 
distributed throughout the galactic volume. How- 
ever, the central concentration of Ha emission in 
ellipticals (e.g. Macchctto et al. 1996) suggests 
that the cooling is concentrated toward the galac- 
tic center where the interstellar density is highest 
and the bulk of the X-ray energy is emitted. 

The motivation of this paper is to explore a 
variety of options for the mass dropout profile of 
cooled gas in bright ellipticals appropriately con- 
strained by the known radial distributions of total 
stellar and non-baryonic mass. The radial mass 
dropout profile of cooled gas cannot be determined 
from first principles because the distribution and 
amplitude of the entropy and magnetic fluctua- 
tions in the hot gas are unknown and difficult to 
evaluate from simple physical arguments. Nev- 
ertheless, the total mass of cooled gas inside an 
effective (half-light) radius Tg must be consistent 
with the mass to light ratio determined from stel- 
lar velocities and with the total mass inferred from 
X-ray observations within re- Assuming that the 
stellar mass to light ratio is uniform with radius, 
the stellar mass M* (r) and the X-ray mass M^ (r) 
appear to be in nearly perfect agreement for two 
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bright Virgo ellipticals in the range O.lrg ^ r < re 
(Brighenti & Mathews 1997a). Because of con- 
straints on the mass distribution of cooled gas pro- 
vided by X-ray and stellar dynamical observations, 
galactic cooling flows provide a critical venue for 
testing the physics of mass deposition in cooling 
flows. 

The mass of cold (T <> 10'' K) gas M^g actu- 
ally observed in ellipticals is many orders of mag- 
nitude less than the total cooled mass estimated 
above. For example, neither HI nor H2 gas has 
been observed in NGC 4472, only upper limits, 
Meg < 10"^ Mq (Bregman, Roberts & Giovanelh 
1988; Braine, Henkel & Wiklind 1988). If stars 
form they must either be indistinguishable from 
the old stellar population or non-luminous. We 
explore here the possibility that most of the cooled 
gas forms dense baryonic clouds or stars that are 
dark at optical and radio frequencies. The very 
low mass stars advanced by Ferland, Fabian & 
Johnstone (1994) satisfy this invisibility criterion, 
while the star formation models of Mathews & 
Brighenti (1999) indicate that (luminous) stars of 
mass ~ 1 — 2 Mq can form in galactic cooling 
flows. 

X-ray studies indicate significant masses of 
cold, absorbing gas in cluster and galactic cool- 
ing flows (e.g. White et al. 1991; Allen et al. 
1993; Fabian et al. 1994; Allen & Fabian 1997; 
Buote 1999), but these results are inconsistent 
with the absence of radio frequency emission from 
the cold gas (Braine & Dupraz 1994; Donahue 
& Vogt 1997) and should be regarded as contro- 
versial until this inconsistency is resolved. Even 
taken at face value, the total mass of cooled gas 
implied by the X-ray absorption in cluster cooling 
flows is typically only a small fraction of the total 
mass that should have cooled in a Hubble time 
(Allen & Fabian 1997: Wise & Saraziu 1999), 
implying that most of the cooled gas may have 
formed stars. If cooled gas forms into small stars, 
these stars will have apogalatica near their point 
of origin where they will spend most of their or- 
bital time. We shall assume that the gas mass 
that cools and drops out of the flow contributes 
optically dark (stellar) mass at the radius where 
the cooling occurred. 

In the following we describe a series of gas dy- 
namical calculations for the evolution of X-ray 
emitting interstellar gas over the Hubble time and 



investigate a variety of assumptions about the ra- 
dial distribution of optically dark cooled gas. To 
be speciflc, we compare our models with the well- 
observed massive elliptical NGC 4472. We find 
that the mass of cooled gas contributes signifi- 
cantly to dynamical mass to light determinations 
within Ve based on stellar velocities. 

2. KNOWN STELLAR AND DARK MASS 
DISTRIBUTION IN NGC 4472 

The E2 elliptical NGC 4472 is a luminous, 
slowly rotating [(i>/o-)* = 0.43] galaxy in the Virgo 
cluster. With an adopted distance d = 17 Mpc, 
its optical luminosity is L_b = 7.89 x 10^° Lbq 
and its half light or effective radius re = 1.733' 
is 8.57 kpc (Faber et al. 1989). The total stel- 
lar mass M^t = 7.26 x 10^^ M© is found from 
the mass to light ratio M/Lb = 9.2 determined 
by van der Marel (1991) with a two-integral stel- 
lar distribution function. This mass to light ratio 
is appropriate to the galactic region within about 
0.4re where stellar velocities are well determined, 
although the mass determined from X-ray obser- 
vations suggests that M/Lb remains constant to 
at least re (Brighenti &: Mathews 1997a). UM/Lb 
is spatially constant, the stellar mass also has a de 
Vaucouleurs profile. Within a central core or break 
radius rj, = 2.41" (200 pc) the stellar density pro- 
file fiattens (Faber et al. 1997), but we shall not 
consider this small feature here. NGC 4472 con- 
tains a central black hole of mass M^h = 2.9 x 10^ 
Mq (Magorrian et al. 1998). 

The total mass distribution in luminous ellip- 
ticals can most easily be determined from the ra- 
dial variation of density and temperature in the 
hot interstellar gas, assuming hydrostatic equilib- 
rium. Figure 1 illustrates the interstellar density 
and temperature profiles in NGC 4472. The filled 
circles in Figure 1 are Einstein HRI observations 
(Trinchicri, Fabbiano & Canizares 1986) and open 
circles are ROSAT HRI and PSPC data from Irwin 
& Sarazin (1996) (also see Forman et al. 1993). 
The T(r) and n{r) profiles have been fit with ana- 
lytic curves as described by Brighenti & Mathews 
(1997a). 

Hydrostatic equilibrium in the hot interstellar 
gas is an excellent approximation since the cool- 
ing flow velocity is very subsonic. The total mass 
interior to radius r determined from X-ray obser- 
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vations is 

M (r) = -^Iil}L f^}^ + t^logT ^ ^ dlogPm \ 
^ G^mp \dlogr dlogr P dlogr / 

(1) 

where nip is the proton mass and fi — 0.61 
is the mean molecular weight. The last term, 
representing the possibility of magnetic pressure 
Pm = B^/Stt, is negative if dPm/dr < as seems 
likely. If the magnetic term is important but not 
included, the total mass will be underestimated. 
Assuming P„i = 0, the total mass Mx{r) in NGC 
4472 is shown with a solid line in Figure lb. 

In the outer halo, r > re, the total mass is domi- 
nated by the dark halo. The dark halo mass distri- 
bution in NGC 4472 can be approximated with an 
NFW halo profile (Navarro, Frcnk & White 1996) 
of virial mass Mh = 4 x 10^^ Mq although the 
observed halo is somewhat less centrally peaked 
than NFW (Brighenti & Mathews 1999aj. Within 
re the contribution of the dark NFW halo mass in 
the model is small; for example at r < r^/S the 
total mass to light ratio is M/ Lb = 10.18, only 10 
percent greater than the dynamic value 9.2 deter- 
mined in r < 0.4re. 

It is remarkable that the total mass found from 
the X-ray data Mx{r) is nearly identical to the 
expected dynamical mass M^{r) (based on stellar 
velocities and M/L) in the range O.lre ^ r ^ re 
(Figure 1). An almost identical agreement in 
this radius range is indicated by X-ray observa- 
tions of another bright Virgo elliptical, NGC 4649 
(Brighenti & Mathews 1997a). In this important 
region the hot gas is confined by the stellar po- 
tential. The excellent agreement of the stellar and 
X-ray masses supports the consistency of two rad- 
ically different mass determinations: from stellar 
velocities and from the radial equilibrium of hot 
interstellar gas. The apparent agreement of the 
X-ray and stellar masses in this range of galactic 
radii also indicates that the hydrostatic support 
of the hot gas is not strongly influenced by local 
magnetic fields and rotation. 

However, it is not obvious why M/Lb would be 
constant with galactic radius, particularly when 
the cooling dropout mass is considered, and why 
the agreement between stellar dynamic and X-ray 
masses no longer obtains in the central regions 
r ^ O.lrg. In this central region the total mass 
indicated by the X-ray observations in Figure 1 



is considerably less than the expected mass based 
on an assumed de Vaucouleurs profile and con- 
stant mass to light ratio. This type of deviation 
could be due to magnetic or other non-thermal 
pressure in this region, to rotation, or to local cool- 
ing dropout in the hot interstellar gas. The lower 
mean temperature in cooling regions lowers the 
total apparent gas temperature and results in an 
underestimate of the total internal mass (equation 
l). 

We assume that currently available Einstein 
and ROSAT observations are accurate in the cen- 
tral region of NGC 4472, r < O.lre. These ob- 
servations have been reduced assuming no (non- 
Galactic) photoelectric absorption by low temper- 
ature gas in the central regions. If X-ray absorp- 
tion is present, the true hot gas density would 
be more centrally peaked than shown in Figure 
1 and the total mass indicated by the X-ray ob- 
servations would increase. Buote (1999) finds that 
two-temperature models fit the X-ray spectrum for 
NGC 4472 quite well. The two temperatures do 
not necessarily need to be spatially mixed; they 
could also approximately represent the range of 
the radial temperature variation observed in NGC 
4472. In Buote's two temperature model, only the 
cooler component (T ~ 0.7 keV located in r ;^ re) 
requires an absorption column Nh = 2.9 x 10^^ 
cm~^ in excess of the Galactic value. However, the 
influence of cold gas having this column density 
on the hot gas density plotted in Figure 1 is small. 
For a worst case example, suppose that absorb- 
ing material with column density Nh = 2.9 x 10^^ 
cm~^ is in a disk oriented perpendicular to the 
line of sight and that this disk absorbs all X-rays 
from the back side of the galaxy. The radius of 
this disk would be < 370 pc if it contained the 
maximum mass Meg ~ 10'' Mq allowed by CO 
and HI observations (Bregman, Roberts & Gio- 
vaneUi 1988; Braine, Henkel & Wikhnd 1997) or 
20 kpc if it contained all of the gas that has cooled, 
Meg = 3 X 10^° Mq. The X-ray surface brightness 
within the opaque disk would be reduced by 2, but 
the corresponding gas density would be lowered 
by only 2^/2 = iff. in 

since the volume emissivity 
oc ri^ . Such a small correction in the density (gra- 
dient) could not account for the large mass dis- 
crepancy between the X-ray mass and the stellar 
dynamical mass shown in Figure 1 within O.lre. 
The densities and temperatures in Figure 1 
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were determined from X-ray data assuming the 
abundance of the hot gas is uniformly solar. 
Since the gas is likely to be more metal rich at 
smaller galactic radii (Matsushita 1997; Brighcnti 
& Mathews 1999b), an allowance for this gradient 
would tend to lower the derived density gradient 
and the internal mass, increasing the discrepancy 
in r ^ 0.1 Te in Figure 1 by a small amount. In the 
following discussion we shall ignore the relatively 
small possible influence of absorption or mctallic- 
ity gradients on the results shown in Figure 1. 

3. HYDRODYNAMICAL MODELS 

The hydrodynamical models we use in this pa- 
per are similar to those in our recent papers (e.g. 
Brighenti & Mathews 1999a) so we provide only a 
brief review here. Hot interstellar gas in ellipticals 
has a dual origin: (i) mass loss from an evolv- 
ing old stellar population and (ii) secondary infall 
into the overdensity perturbation that formed the 
galaxy group within which the elliptical formed 
by early merging events. For a given set of cos- 
mological parameters, dark and baryonic matter 
flow toward an overdensity region. The dark mat- 
ter forms an NFW halo, growing in size with 
time. Spherical geometry is assumed. Within 
the accretion shock at time = 2 Gyr, when 
enough baryons have accumulated, we form the 
current de Vaucouleurs stellar profile and release 
the energy of all Type II supernovae according 
to a Salpeter IMF (slope: x = 1.35, mass lim- 
its: nn = 0.08 and m„ = 100 Mq). All stars 
greater than 8 Mq produce Type II supernovae 
each of energy Egn = 10^^ ergs. We assume that 
a fraction egn = 0.8 of this energy is delivered to 
the internal energy of the gas. We have shown 
(Brighenti & Mathews 1999b) that such a galaxy 
formation scheme can work in a variety of cos- 
mologies: flat or low density, with or without a 
lambda term. The evolution of gas within the op- 
tical effective radius r^, of most interest here, is 
insensitive to these cosmological parameters. For 
simplicity therefore we assume a simple flat uni- 
verse with = 1, Ho = 50 km s~^ Mpc~^ and 
fif, = 0.05. We characterize the dark halo with 
an NFW proflle (Navarro, Frenk & White 1996) 
having a virial mass Mh = 4 x 10"'^^ Mq at the 
current time t„ = 13 Gyrs. The models we dis- 
cuss here are identical to the standard model of 
(Brighenti & Mathews 1999a) except we now use 



a finer central spatial zoning (65 pc for the inner- 
most zone), a SNII efficiency = 0.8, and a mass 
"dropout" function q{r) with more adjustable pa- 
rameters (see below). 

Our objective is to seek solutions of the gas- 
dynamical equations including mass dropout that 
jointly satisfy several observational constraints at 
time tn- (i) the observed hot gas density, tempera- 
ture and X-ray surface brightness profiles, (ii) the 
known dynamical mass in the galactic center usu- 
ally attributed to a massive black hole, (iii) the ap- 
parent dynamical mass to light ratio M/Lb = 9.2 
determined in r ^ 0.4re, and (iv) an apparent 
internal mass M^ir) in (0.1 — l)re based on Equa- 
tion (1) that agrees with the constant- (M/L^) de 
Vaucouleurs profile as shown in Figure 1. 

The baryonic component in our models has a 
complex evolution. Much of the initial baryonic 
mass is consumed in creating the stellar system. 
When the Type II supernova energy is released, 
a significant mass of gas is expelled as a galac- 
tic wind. After these early events, the interstel- 
lar gas is re-established and sustained by stellar 
mass loss and by inflow of circumgalactic gas (sec- 
ondary infall), most of which was previously en- 
riched and expelled by SNII. We assume that the 
stars form during a short epoch that can be de- 
scribed by a single burst Salpeter IMF as dis- 
cussed above. The stellar mass loss rate for this 
IMF varies as = Q!*(t)M*t where a*(i) = 

4.7xl0-2O[t/(f„-t^^)]-i-3 sec-^ Although galac- 
tic stars form at i*^ = 1 Gyr, their mass loss con- 
tribution to the ambient interstellar gas is assumed 
to begin at a later time, = 2 Gyrs. Since galac- 
tic stars have been enriched by supernova ejecta, 
the single burst model cannot be strictly correct, 
but our approximation Q:*(f < t*) = is consis- 
tent with several early starbursts c;loscly spaced in 
time and allows for metal enrichment of old galac- 
tic stars that are not in the flrst single-burst pop- 
ulation. If the de Vaucouleurs profile is a result 
of largely dissipationless merging, some or most of 
the star formation must have occurred at a time 
t<,s before the important merging events at = 2 
Gyrs. By taking t^s < t* we reduce by ~ 10^° Mq 
the total amount of gas ejected by stars within 
the galactic potential. We recognize the inconsis- 
tencies in these approximations of complex stellar 
formation and dynamical processes that are poorly 
understood. However, once the galaxy is formed. 
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we follow the interstellar gas dynamical evolution 
in full detail, conserving mass and energy. 

Continued heating by Type la supernova is as- 
sumed to vary inversely with time, SNu(i) = 
SNu(t„)(t„/t), where the current rate, SNu(t„) = 
0.03 SNIa per 100 yrs per 10^° is©, is near the 
lower limit of observed values, SNu(i„) = 0.06 ± 
0.03(iJ/50)2 (Cappellaro et al. 1997), as required 
to maintain the low interstellar iron abundance. 

For the models discussed here the equation of 
continuity includes a "mass dropout" term: 



sider a wider range of constant or variable dropout 
coefRcients q{r) parameterized by 



dp 
dt 



or ^ tdo 



where tdo = 5mpkT/2iJ,pA is the time for gas to 
cool locally by radiative losses at constant pressure 
(see, c. g. Sarazin & Ashe 1989). The cooling is 
assumed to be instantaneous without advection in 
the cooling flow, i.e. tdo ^ tfiow = r/v, although 
in practice this inequality may not always be sat- 
isfied. While this type of cooling dropout has been 
widely used in past models, there is no adequate 
physical model for mass dropout. Clearly, the gas 
must cool somewhere - the emission of X-rays in- 
dicates a large net energy (and mass) loss from the 
interstellar medium. 

When small regions of low entropy cool, the 
pressure remains constant since the sound cross- 
ing time is much less than the flow time. Fol- 
lowing Fabian, Nulsen & Canizares (1982) and 
Ferland, Fabian, & Johnstone (1994), we assume 
that cooled gas converts to a second population 
of optically dark, low mass stars. Regarding Ha 
emission as a tracer for the cooling gas, Math- 
ews & Brighenti (1999) have shown that the cool- 
ing occurs at a multitude (~ 10^) of cooling sites 
distributed throughout the inner galaxy and that 
only stars with mass ;^ 1 — 2M© can form at each 
cooling site. This is supported by the observed 
absence of young massive stars and SNII in ellip- 
tical galaxies. Nevertheless, in the following dis- 
cussion we assume that the maximum stellar mass 
in the dropout population is sufficiently low that 
the optical light from these stars is unobservable. 
We expect cooling and low mass star formation 
to be concentrated within at least 2 kpc, which is 
the observed extent of Ha emission in NGC 4472 
(Macchetto et al. 1996). While the region con- 
taining optical emission lines provides a natural 
guideline for selecting the dropout profile, we con- 



q{r) = qoex^^i-r/rdoY 



(2) 



which concentrates the cooling within radius Tdo- 
Note that even when q is constant, the mass 
dropout term is spatially concentrated, p/tdo oc 



P'. 



4. MODELS WITH ZERO OR CON- 
STANT q 

In a series of recent papers in this Journal 

we have presented evolutionary cooling flow mod- 
els for NGC 4472 that agree quite well with the 
observed distributions of interstellar gas density, 
temperatiirc and metallicity (Brighenti & Math- 
ews 1999a; 1999b), particularly at intermediate 
and large radii. Since the total mass of cooled 
gas in these calculations was only a few percent of 
the stellar mass, the gravitational contribution of 
cooled gas was ignored, even in models including 
mass dropout. For a better understanding of the 
inner galaxy, we now include the gravity of all gas, 
both hot and cold, except when specifically noted. 
In this section we begin with models having q = 
everywhere, so that cooling to low temperatures 
occurs in a small central region, then we investi- 
gate models in which q is constant throughout the 
cooling flow. 

4.1. Models Without Distributed Mass Dropout 

(q = o) 

We begin by considering a perfectly homoge- 
neous interstellar flow in which all the gas reaches 
the central computational zone (r = 65 pc) or its 
neighboring zones where it then cools to T <C lO'' 
K. For comparison we discuss two cases: in Model 

1 we ignore the gravity of this cooled gas; in Model 

2 (and subsequent models) we include its gravi- 
tational influence on the flow. The mass of gas 
that cools into the central gridzone is dynami- 
cally equivalent to a massive black hole, but we 
do not necessarily regard our calculation as a real- 
istic model for black hole formation. In Figure 2a 
we illustrate the radial variation of gas density and 
temperature for Models 1 and 2 after the interstel- 
lar gas has evolved to t = i„ = 13 Gyrs. Several 
global parameters for these and subsequent mod- 
els are listed in Table 1. 
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When the gravity of the cooled gas is considered 
(Model 2), the interstellar gas within a few kpc 
of the galactic center is compressed and sharply 
heated in the local potential. Such hot central 
thermal cores are not generally observed (but see 
Colbert & Mushotzky 1998). The central mass in 
both models, Mcent = 4.65 x 10^" Mq is about 
16 times larger than the black hole mass found in 
NGC 4472, Mbh = 2.9 x lO'' Mq (Magorrian et 
al. 1998). For those reasons neither Model 1 nor 
2 provides a realistic interstellar mass distribution 
for this galaxy. 

ROSAT band X-ray surface brightness profiles 
T,x{R) for Models 1 and 2 are shown in Figure 
2b and the total ROSAT band luminosities are 
listed in Table 1. For Model 1 the X-ray bright- 
ness peaks strongly in the galactic center, diverg- 
ing from the observations within about 3 kpc; it 
was just this sort of disagreement that initially 
led to the assumption of distributed mass dropout 
in cooling flows (Thomas 1986; Vedder, Trester, 
& Canizaros 1988; Sarazin & Ashe 1989). If the 
efficiency of heating by early SNII were lowered. 
Models 1 and 2 could be made to agree better with 
T,x observations in r > 10 kpc, but the computed 
would rise even further above the observations 
within a few kpc. The X-ray luminosities for 
Models 1 and 2 listed in Table 1 are unreliable in 
part because of numerical inaccuracies caused by 
the extremely steep variation of gas temperature 
and density in the central 2 or 3 computational 
zones. This numerical difficulty is common to all 
cooling flows that proceed to the very center of 
the galaxy before cooling, for any reasonable cen- 
tral grid spacing. Since more Pdv work is done on 
the flow in Model 2 (where the gravity of cooled 
gas is included) , we expect that should also be 
larger than for Model 1. The opposite sense of 
the change in shown in Table 1 evidently re- 
sults from numerical inaccuracies near the central 
singularity where zone to zone variations are no 
longer linear. If the gas has not cooled before flow- 
ing into the core (< 100 pc), as in Models 1 and 
2, a significant fraction of the total X-ray emission 
should come from this very central region, again 
in disagreement with observations. 

In Models 1 and 2 a sphere of cold (T = 10 K), 
dense gas accumulates in the central gridzones and 
grows in mass and size over the Hubble time. This 
sphere is an unrealistic artifact of our computa- 



tional assumptions. Therefore, to explore further 
the central numerical difficulty with encoun- 
tered in Models 1 and 2, we considered two addi- 
tional models using higher spatial resolution (ra- 
dius of central zone is only 15 pc), both of which 
include the gravitational influence of cooled gas. 
The first model (Model 2.1) is similar to Model 
2 with g = in all zones. In the second model 
(Model 2.2) we set (/ = in all zones except the 
central zone where q = 4. Model 2.2 is the ap- 
propriate limit of the series of models discussed 
below in which the cooling dropout is more and 
more centrally concentrated. 

While the flow in r > 1 kpc is very similar in 
Models 2.1 and 2.2, the behaviour at smaller radii 
is quite different. For example the flow parame- 
ters for Model 2.1 at 100 pc and tn = 13 Gyrs are: 
T = 7.9 X 10'^ K, n = 1 cm'^, and u = -320 km 
s~^. At the same radius and time for Model 2.2 
the flow is quite different: T = 6.3 x 10"" K, n = 5 
cm~^, and u = —60 km s~^. At projected radius 
R = 100 pc, the ROSAT X-ray surface brightness 
in Model 2.2 is 20 times larger than Model 2.1 
and the total ROSAT X-ray luminosity integrated 
over the entire cooling flow in Model 2.2 is larger 
by a factor 32. The differences between these two 
models result from upstream propagation of in- 
formation about the central boundary conditions 
made possible by subsonic flow near the origin. 
Fortunately, these numerical and physical difficul- 
ties near the origin do not arise in more realistic 
cooling flows discussed below in which g > at 
larger galactic radii. 

4.2. True and Apparent Gas Density and Tem- 
perature 

While solutions of the gasdynamical equations 

provide the temperature as a function of physi- 
cal radius T(r), the observed temperature T{R) is 
an emission- weighted mean temperature along the 
line of sight at projected radius R. For symmetric 
galaxies and at small galactic radii, these two tem- 
peratures are nearly identical because of the steep 
radial variation in X-ray emissivity. The variation 
of temperature with physical radius T(r, tn) in the 
background flow is shown with light solid lines in 
Figure 2a. 

In the presence of spatially distributed cool- 
ing dropout, the local temperature is an emission- 
weighted mean of the background (uncooled) gas 
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and the cooling regions. Cooling is assumed to oc- 
cur at a large number of cooling sites where the gas 
remains in pressure equilibrium as it cools, appar- 
ently imrcstrictcd by magnetic stresses (sec Math- 
ews & Brighenti 1999 for details). The heavy solid 
lines in Figure 2a show the mean apparent tem- 
perature Teff{r,tn) including contributions from 
locally cooling regions, 



T + qTAi(T) 

i + gAo(T) 



(3) 



where T is the background flow temperature and 
the slowly- varying functions Ai{T) are plotted in 
Brighenti & Mathews (1998). Note that the ef- 
fective temperature is independent of the local 
gas density. The temperature that is actually ob- 
served is an average of Teff along the line of sight; 
this temperature Teff{R, tn) is shown with dotted 
lines in Figure 2a. 

Similarly, the apparent hot gas density is 
increased because of additional emission from 
denser, cooling-out gas. The observed electron 
density shown in the Figure 2a is found by Abel 
inversion of the X-ray surface brightness distri- 
bution. When cooling sites are present, the local 
emissivity into the ROSAT energy band is 

£AE = {pefs/mpfK^E{T) 

= {p/mpfhAE{T)[l + gAo(T)] ergs/sec cm^ 
The effective density is therefore 

ne//=n[l + <zAo(r)]i/2 (4) 

where n is the electron density of the back- 
ground, uncooled gas. The observed (azimuthally- 
averaged) densities in NGC 4472 plotted in Fig- 
ure 2a should be compared with ne//(r, tn) shown 
with heavy solid lines. 

4.3. Distributed Mass Dropout with Constant q 

Lacking an acceptable physical model for spa- 
tially distributed mass dropout in galactic cooling 
flows, we are at liberty to choose any variation 
for the dropout coefficient g(r), appropriately con- 
strained a posteriori by the known dynamical mass 
from stellar velocities, the X-ray mass, the central 
black hole mass and the observed radial variation 
of hot gas density, temperature and X-ray surface 
brightness. It is natural to begin with constant 



q solutions similar to those considered by Sarazin 
& Ashe (1989) in their steady state cooling flow 
solutions. For simplicity we assume that the mass 
of cooled gas remains at the cooling site where it 
contributes to the gravitational potential. 

In the central panels of Figure 2a we illustrate 
the hot interstellar gas density and temperature 
that results after t„ = 13 Gyrs assuming uniform 
q = 1 and 4; these are listed in Table 1 as Mod- 
els 3 and 4 respectively. When q is a constant 
independent of radius, the ratios of true to ap- 
parent values - n/ueff, T/Te/j and Mtot/M^ - 
ilso approximately uniform with galactic ra- 
dius. Here Mtot{r) is the true total mass within 
r and Mx{r) is the value that would be deter- 
mined from X-ray observations of the models by 
assuming hydrostatic equilibrium (as in Fig. 1). 
The influence of constant-g' mass dropout is simi- 
lar at all galactic radii since the factors that con- 
vert background temperature and density to effec- 
tive values in equations (3) and (4) depend only 
on T(r) which is slowly varying, not on n(r). Be- 
cause of the enormous volume and time available 
to gas dropping out at large r, the total mass of 
cooled gas (listed in Table 1) becomes very large 
in the outer galaxy. In Figure 3 we compare the 
radial distribution of stellar and dark halo mass 
with the distributed dropout mass that results af- 
ter 13 Gyrs with g = 1 and 4. The de Vaucouleurs 
"stellar" mass profile in Figure 3 is constructed 
assuming uniform M/Lb = 9.2. 

As q increases from 1 to 4, the background hot 
gas density (light solid lines in Fig. 2a) decreases 
and its radial gradient flattens, causing a rise in 
temperature to provide enough prcssTire to sup- 
port the cooling flow atmosphere. Note that the 
total dropped-out mass within re/3 decreases with 
larger q (Figure 3 and Table 1). Although more 
overall cooling dropout occurs as q increases, most 
of this cooling occurs at very large r and, some- 
what paradoxically, less gas remains for dropout 
closer to the galactic center, r < r^. However, the 
effective (i.e. apparent) density (heavy solid lines 
in Fig. 2a) is less sensitive to q. The q — 1 solu- 
tion (Model 3) is preferred for its fit to observed 
tempcratiires while the q = A solution (Model 4) 
agrees better with the observed density in r < 10 
kpc and almost exactly with the X-ray surface 
brightness (Fig. 2b). 

In addition, the centrally concentrated dropped 
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out mass in the q — 1 solution (Model 3) con- 
tributes a larger fraction of the dynamical mass 
in r ^ re/3 and the total mass within the central 
gridzonc (r = 65 pc) is almost equal to the known 
mass of the black hole in NGC 4472 (Table 1). Ev- 
idently, cooled masses in models with q <1 would 
exceed the central dark mass observed. 

For both values of q considered, the dropped 
out mass listed in Table 1 contributes appreciably 
to the total mass within r = re/3. When both old 
stellar and dropout mass are included, the M/Lq 
values at re/3 shown in Table 1 - 12.19 and 11.21 
- exceed both the dynamical value M/Lb = 9.2 
found by van der Marel (1991) within ^ 0.4re 
and the value M/LB(re/3) = 10.18 in our mod- 
els which includes a small additional contribution 
of non-baryonic dark matter. To quantify the in- 
fluence of the dark baryonic dropout mass on the 
mass to light ratio we include in Table 1 an entry 
for 



vK^e/3) 



{M/LB)..dh.do{^e/i) 

{M/LB).Mre/i) 



- 1 



_ M,,dfe,do(re/3) ^ 
M^,dh{re/2>) 

where M*,d/i(re/3) = 1.35 x 10^^ M© is the com- 
bined mass of luminous stars and dark halo mat- 
ter at re/3 and M^^dh,do{re/^) = Mtot{re/^) also 
includes the dropout mass within this radius. 
Since our computed total mass M^^dh,do{re/'i) ex- 
ceeds the observed dynamical value 9.2 (and also 
10.18), to be fully consistent we should have cho- 
sen M/Lb^ < 9.2 for the old stellar population 
(see below), provided M/Lb for stars produced 
from the cooled gas is infinite as we have assumed. 

If these constant q models are physically appro- 
priate, the agreement between Mx{r) and M, (r) 
in (0.1 — l)re, as shown in Figure 1 (and for 
NGC 4649), is surprising since Mx{re/3) and 
M^^dh{re/S) differ by 10 - 40 percent (Table 1) 
for the constant q models. This suggests that 
the mass to light ratio of the dropout stellar pop- 
ulation is not infinite as we have assumed, but 
comparable with that of the old stellar popula- 
tion. 



5. MODELS WITH CENTRALLY CON- 
CENTRATED DROPOUT 

We now seek evolutionary gasdynamical solu- 
tions with variable dropout coefficients q{r) that 
strongly concentrate the mass dropout in the inner 
galaxy, r ^r^. The limited spatial extent of opti- 
cal emission lines in the cores of bright ellipticals 
suggests that cooling dropout occurs well within 
Te. For example, the Ha -I- [Nil] image of NGC 
4472 is observed out to approximately ~ 0.25re or 
2 kpc (Macchetto et al. 1996). As hot interstellar 
gas cools in this region, its temperature pauses at 
T ~ 10^ K where the gas is heated and ionized 
by stellar UV radiation. Therefore, we consider 
parameters rdo and m in equation (2) that con- 
centrate the mass dropout in the central region, 
but not at the very center as in Model 2. 

The dropout parameters ^o, ^do and m.. arc 
listed in Table 1 for Models 5, 6 and 7. The 
mass dropout in these three models is progres- 
sively more concentrated toward the galactic cen- 
ter. The results for Model 5, for which the dropout 
scale length is very large, Vdo = 15 kpc, are simi- 
lar in most respects to those of Model 4 except the 
massive dropout in the outer galaxy is no longer 
present. In particular, the mass to light ratio at 
re/3 in Table 1 is very similar for Models 4 and 5. 

The current interstellar density and tempera- 
ture variations for Model 6 with rdo = 2 kpc are 
shown in Figure 2a and Tix{R) is plotted in Fig- 
ure 2b. The apparent density and Y^x (heavy solid 
lines) are considerably greater than observed val- 
ues in r < 3 kpc, although the projected apparent 
temperature is a reasonable fit to the NGC 4472 
data. The radial distribution of dropout mass for 
Model 6 is shown in Figure 3. From Table 1 the 
total (old stars, halo and dropout) mass to light ra- 
tio at re/3 is Mtot/LB = 13.03, 29 percent higher 
than van der Marcl's value and 22 percent greater 
than the corresponding value for our background 
model galaxy. 

Also shown in Figures 2a and 2b are similar re- 
sults for Model 7 in which the mass dropout, now 
approximated with a Gaussian, is concentrated 
within Vdo = 800 pc. Model 7 is constructed so 
that most of the mass dropout occurs in r ^ O.lre, 
corresponding to the region of apparent disagree- 
ment between and M^, in Figure 1. It is of 
interest that the gas density and S^, within 1 kpc 
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for Model 7, where the dropout is greatest, exceeds 
those of Model 2 in which there is no dropout at 
all. This can be understood from the flow veloc- 
ity distribution. For Model 2 without distributed 
dropout the inward moving gas velocity increases 
through the entire region illustrated and shocks at 
r 1 kpc. In the presence of dropout, the flow 
velocity at corresponding radii in Model 7 is much 
slower and reaches a maximum (negative) value 
near r = 1 kpc then approaches zero subsonically 
at the origin. The density and T,^ enhancements 
in the background flow for Model 7 at r ^ 3 kpc 
are due to a local compression as the gas flows 
into slowly moving gas in the core. Within about 
3 kpc, both the background and apparent densi- 
ties exceed the observations by a larger factor than 
those of Model 6 and T,^ also peaks unrealistically 
in this same region (Fig. 2b). The dropout mass 
for Model 7 shown in Figure 3 equals that of the 
old population stellar mass at r « 1 kpc (0.12re). 
For these reasons Model 7 seems less satisfactory 
than Model 6, but neither is as generally successful 
as Model 3 (g = 1). 

Although increased mass dropout can decrease 
the X-ray surface brightness 'Sx, as when uniform 
q increased from 1 to 4 (Models 3 and 4 in Fig. 2b), 
this is not always the case. When the dropout is 
concentrated more toward the galactic center, T,x 
actually increases, as in the transition from Model 
6 to 7. 

It is interesting to determine the influence of 

distributed dropout on the total apparent mass 
Mx{r) found from the model by assuming hydro- 
static equilibrium (equation 1). Due to the contri- 
bution of low temperature cooling regions to the 
total X-ray emission, T{r) and therefore Mj:(r) is 
always lower than the true mass M^_dh{r) in cool- 
ing dropout regions. A difference in this sense is 
apparent in Table 1 at r = re/3 for Models 3-7; 
this is similar to the mass discrepancy in Figure 1 
at r ^ O.lrg. 

The high central apparent gas density and sur- 
face brightness for Model 6 shown in Figures 2a 
and 2b are obvious problems for this model. How- 
ever, the gas density can be reduced if a strong 
magnetic fleld or other non-thermal energy den- 
sity is present in r < 0.25re = 2 kpc. Dynamically 
important magnetic flelds may also be required 
to fit the X-ray data of NGC 4636 (Brighenti & 
Mathews 1997a) and may be generally expected in 



luminous ellipticals (Mathews & Brighenti 1997; 
Godon, Soker & White 1998). Additional non- 
thermal pressure support is also implied for Model 
6 since the effect of central cooling dropout fails 
to lower the apparent mass Mx{r) below the ac- 
tual mass Mtot as much as the observed deviation 
shown in Figure 1. Like many bright ellipticals, 
NGC 4472 has a weak non-thermal radio source 
within the central ~ 4 kpc, indicating B ~ 10—100 
fiG (Ekers & Kolanyi 1978). 

For all calculated models in Table 1 with dis- 
tributed dropout, the total mass Mtot = Mt,^dh,do 
significantly exceeds the mass of the old stel- 
lar population plus dark halo M^^^h throughout 
the inner galaxy, indicating that dropout material 
makes an important additional contribution to the 
total mass. However, the apparent mass M^, de- 
termined from equation (1) is less than Mtot in 
the inner galaxy and, for Model 6, can also be 
less than M*. Of particular interest is the region 
O.lre ^ r < re (i.e., -0.07 logrupc ^ 0.93) 
in Figure 1. Although the agreement in Figure 
1 is excellent in this region, for Models 6 and 7 
the total mass is larger than M^,^d/j(r) and values 
of A„/; in Table 1 suggest that the dropped out 
mass contributes 25 - 35 percent of the total mass 
in this region. Therefore, if the dropout mass is 
optically dark, the true stellar mass to light ratio 
of the old stellar population must be M/Lb « 6 
rather than the value M/Lb = 9.2 found by van 
der Marel which includes the dropout mass. 

To investigate such a possibly more self- 
consistent old stellar component, we consider 
Model 8 based on the same q{r) used in Model 
6, but with M/Lb = 6 for the old stellar pop- 
ulation. For additional consistency in Model 8, 
a*(t) is increased by the ratio of assumed stellar 
mass to light ratios 9.2/6 = 1.53 as described be- 
low; the total stellar mass ejected is identical to 
that in Model 6. Figure 3 illustrates the dropout 
mass profile for Model 8. The apparent density, 
temperature and S^, profiles for Model 8 shown 
in Figures 4a and 4b are very similar to those of 
Model 6, so this adjustment of the stellar M/Lb 
has had little effect. 

The radial mass profiles of Models 6 and 8 
are compared in Figure 5. In this plot the open 
circles show the observed X-ray mass Mx{r) for 
NGC 4472 using equation (1) and the solid lines 
shows Mx{r) based on equation (1) using n(r) and 
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T(r) from the models. The dashed hnes show 
Mtot{r) and the dotted hnes are the stellar mass 
M*(r) based on a de Vaucouleurs profile with 
M/Lb = 9.2 in the upper panel (Model 6) and 
M/Lb = 6 in the lower panel (Model 8). The 
superiority of Model 8 is evident from the closer 
agreement between the X-ray mass data points for 
NGC 4472 and the solid line for that model. This 
agreement for Model 8 would be even closer in 
the range logr « 0.5 — 1 if we had used a dark 
halo mass profile less centrally peaked than NFW. 
Model 8 may provide the most self-consistent over- 
all fit to the mass constraints for NGC 4472; if so, 
the old stars in NGC 4472 have a mass to light 
ratio M/Lb ~ 6, about ^ 30 percent lower than 
the mass to light ratio determined from stellar dy- 
namics. 

In summary, none of our models is fully satis- 
factory in every respect. The total mass to light 
ratio (including dropout mass) in r ^ re/3 is 10 
- 35 percent higher than the value for the un- 
derlying galaxy. However, in Model 8 in which 
M/Lb = 6 for the old stars, the difference be- 
tween the X-ray mass determined from the mod- 
els and NGC 4472 is appreciably reduced. Nev- 
ertheless, the central apparent gas density and X- 
ray surface brightness in Model 8 are still larger 
than observed, requiring additional non-thermal 
support. There is no independent theoretical jus- 
tification for the dropout profile q{r) assumed in 
Models 6 and 8; as explained earlier, the dropout 
distribution depends on unknown interstellar en- 
tropy fluctuations. Nevertheless, for all models 
considered here the mass of cooled interstellar gas 
contributes significantly to the total mass and the 
dynamically determined mass to light ratio within 
the inner galaxy. 

6. FINAL REMARKS AND CONCLU- 
SIONS 

In this series of calculations we have taken a 
census of all baryons involved in the evolution of 
a large elliptical galaxy: the original stellar com- 
ponent, the interstellar medium, and - of most in- 
terest - the small but troublesome mass of hot gas 
that cools over cosmic time. We have shown that 
the radial distribution of cooled interstellar gas in- 
fluences dynamical and X-ray determinations of 
the total interior mass and the radial profiles of 



apparent density, temperature and X-ray bright- 
ness of the hot gas. In our models, cooled gas is 
slowly deposited in the central galaxy r ^ rg (see 
Figure 3) as indicated by the extent of observed 
Ha emission. Since there is little or no direct ob- 
servational evidence for the mass that has dropped 
out in ellipticals like NGC 4472, the cooled mass 
must either be dark at optical and radio frequen- 
cies or indistinguishable from the old stars. Low 
mass stars are an obvious and physically reason- 
able endstate for the cooled gas (Ferland, Fabian 
& Johnstone 1994; Mathews & Brighenti 1999). 
We also suppose that the cooled gas remains at 
the dropout site where it contributes to the galac- 
tic potential. 

6.1. Reducing Stellar Mass Loss 

In an attempt to reduce the infiuence of cool- 
ing and cooled gas on the models, we have altered 
many of the model parameters. The chosen cos- 
mology (fi = 1; n = 0.3 and Q\ = 0.7, etc.) 
or the baryon fraction f2(, have little influence on 
the total dropout mass. Changing the times when 
the stars and the galactic potential form (t*s and 
tt.) or the spatial scale of the release of SNII en- 
ergy within reasonable limits have only a modest 
infiuence on the dropout mass that accumulates 
by time t„ = 13 Gyrs. Increasing the interval 
t*s — between star and galaxy formation does 
reduce the total dropout mass, but this interval 
cannot be too large since luminous ellipticals are 
observed at large redshifts. The cooling dropout 
would not be dramatically reduced if massive el- 
lipticals were all only a few Gyrs old since most of 
the stellar mass loss occurs just after star forma- 
tion; however, many or most luminous ellipticals 
are thought to be very old. 

Perhaps the most effective way to preserve the 
excellent agreement between and M^ in Fig- 
ure 1 in O.lre ^ r ^ Ivg, without constraining 
the mass dropout profile q{r), would be to re- 
duce the total mass that has cooled over cosmic 
time. The most sensitive parameter influencing 
the cooled mass is the specific rate of stellar mass 
loss, a*(t). Although the total mass of hot gas 
increases with time due to the continued influx of 
secondary infalling gas, the mass of hot gas within 
the optical galaxy originates mostly from stellar 
mass loss and the X-ray luminosity there scales as 
Lx oc a*(t) oc t~'^-^ (Appendix B of Tsai & Math- 
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ews 1995). Computed models similar to those de- 
scribed here but with arbitrarily reduced a* (t) fit 
the n(r) and T,x{R) data rather well at time tn 
and produce much less cooling dropout. However, 
a*(t) cannot be lowered without also increasing 
the stellar mass. For all reasonable power law ini- 
tial mass functions, a,(t) varies inversely with the 
stellar mass to light ratio: 



Q!*(i) 



dM^/dt 



dM^/dt 



The physical explanation for the constancy of 
{dM^/dt)/ Lb is that both Lb (dominated by 
post-main sequence stars) and dM^/dt depend on 
the instantaneous rate that stars leave the main se- 
quence so this IMF-dependent factor cancels out 
(e.g. Rcnzini & Buzzoni 1986). To illustrate this 
result, we use the Renzini-Buzzoni procedure and 
plot in Figure 6 the relationship between a* and 
M/Lb at time t„ = 13 Gyrs for 84 power law 
IMFs with slopes x =0.6, 0.8, 1.0, 1.2, 1.4, 1.6, 
and 1.8, each with lower and upper mass limits of 
mi =0.01, 0.032, 0.1, 0.32 and m„ =10, 32, 100. 
The remarkable linearity in Figure 6 shows that 
the quantity in square brackets in the equation 
above is almost invariant to large changes in the 
slope or mass limits for power law IMFs. 

Therefore, if a^{t) is reduced by 2 or 3 in an at- 
tempt to reduce the total dropout mass, the stel- 
lar mass (and M/Lb) must be increased by the 
same factor; as a result the total mass ejected, 
/ a^M^^dt, and the total dropout mass do not 
change. 

6.2. Effect of Galactic Rotation 

Although massive ellipticals are not rotation- 
ally flattened, they do rotate significantly, e.g. 
(w/cr)* = 0.43 for NGC 4472 (Faber et al. 1997). 
If the hot interstellar gas rotates in the same sense 
as the bulk of mass-losing stars, stars formed from 
the cooled gas should form into a disk of scale ~ re 
(Brighcnti & Mathews 1997b), although the devel- 
opment of such disks is likely to be suppressed by 
the mass dropout process. Nevertheless, to the 
extent that the cooled gas has a disk-like distribu- 
tion, its global influence on the stellar dynamics in 
r ^ 0.4re would be less than if the same dark mass 
were distributed spherically. Remarkably, there is 
no observational evidence at present for rotational 



flattening in the X-ray images of giant ellipticals 
like NGC 4472. 

6.3. Non-Baryonic Dark Matter within r^ 

Dynamical determinations of the mass to light 
ratio from stellar velocities reflect the entire mass 
within the stellar orbits, including non-baryonic 
mass. The NFW halo we use in our models for 
NGC 4472 agrees with X-ray observations in the 
extended halo but is slightly too massive near 
r ~ Te relative to the X-ray mass Mx{r) obser- 
vationally determined for NGC 4472. This may 
indicate that dark halos are less centrally peaked 
than NFW (see Kravtsov et al. 1998). The mass 
of our NFW halo model contributes about 10 per- 
cent to the dynamical mass to light ratio mea- 
sured within re/3, but the NFW profile is proba- 
bly disturbed in this region. When the dominant 
baryonic mass in r ;^ re compressed to form the 
de Vaucouleurs profile, we expect that the NFW 
halo was dragged inward and distorted. How- 
ever, the dark halo cores of the earlier galactic 
condensations that merged to form the elliptical 
may have expanded due to starburst driven galac- 
tic winds. Because of these various counteracting 
effects, the small non-baryonic contribution to dy- 
namical mass to light determinations is uncertain. 

6.4. Contribution of Dropout to "Stellar" M/L 

For all models studied here - based on a wide 
variety of mass dropout profiles q{r) - the mass 
of cooled interstellar gas contributes substantially 
to the total mass within ^ QAr^ where the stellar 
mass to light ratio is determined from stellar ve- 
locities. If optically dark low mass stars form from 
the cooled gas, the "stellar" mass to light ratios in 
the literature refer to two distinct stellar popula- 
tions having radically different initial mass func- 
tions and spatial distributions. The stellar mass 
to light ratio of the original, optically luminous 
stellar population is lower than published values 
indicate. 

A complete solution of this problem requires a 
better understanding of the physics of star for- 
mation and the processes that control the stel- 
lar IMF. We have assumed here that a Salpeter 
IMF provides a satisfactory approximation to the 
original single-burst star formation at early times. 
Yet we argue that the younger dropout IMF is 
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strongly skewed toward low mass stars. The IMF 
has evolved over time. It is possible therefore that 
the early, more intense mass dropout resulted in a 
more nearly Salpctcr-likc IMF, producing a frac- 
tion of currently observed luminous stars in ellip- 
ticals. If so, this early dropout would not con- 
tribute to the excess dropout mass that wc find 
in our models, but would have already been in- 
cluded in the original de Vaucouleurs population. 
High density, metal enriched stellar cores in ellip- 
tical galaxies may have derived from normal-IMF 
star formation from early, more intense interstellar 
cooling. This is similar to assumptions made for 
dissipative galactic core formation from the con- 
vergence of gas following major mergers (Mihos 
& Hernquist 1996). While such notions cannot be 
entirely dismissed, the approximate universality of 
the de Vaucouleurs light profile among ellipticals 
may argue against a dual formation process for 
the radial distribution of luminous stars: violent 
relaxation and cooling flow dropout. 

Throughout this discussion we have assumed 
that the stellar population formed from cooling 
flow dropout is optically dark. Although there 
is little or no evidence that normal massive OB 
stars (or SNII) are present in elliptical galaxies, it 
is possible that younger stars having masses up to 
^ 1 — 2 Mq are present and that such intermediate 
mass stars could form from the cooled gas (Math- 
ews & Brighenti 1999). This type of dropout stel- 
lar population could contribute to the total optical 
light. If the mass to light ratio of dropout and old 
stellar populations are similar, the dropout com- 
ponent could be difflcult to detect by the means 
we have discussed here and its perturbation on 
the observed M/Lb would be greatly lessened. In 
this case the dropout population would introduce 
an additional radial light profile that would dif- 
fer slightly from that of the old stellar population. 
The dropout mass profiles in Figure 3 indicate that 
Models 4 and 6 would be rather diflicult to detect 
against the background stellar light. Intermediate 
mass dropout stars could therefore provide a sat- 
isfactory resolution to the problems we have dis- 
cussed here. 

6.5. Influence of Dropout on the Fundamental 
Plane 

The ensemble of elliptical galaxies is known to 
have global parameters that deviate slightly from 



the assumptions of virial equilibrium and homol- 
ogous structure. The deviation of this fundamen- 
tal plane relationship is in the sense that the dy- 
namical mass to light ratio increases with galactic 
mass, M/Lb oc M°-^^ (Dressier et al. 1987; Djor- 
govski & Davis 1987). Such a non- homologous 
deviation could in principle be produced by the 
small amount of cooled interstellar gas Meg pro- 
vided it increases appropriately with M*. To test 
this idea, we performed an identical hydrodynamic 
calculation for an elliptical galaxy having a mass 
one fourth that of NGC 4472. The dark halo mass 
was also reduced by the same factor but the cos- 
mological environment and mass dropout distri- 
bution were identical to those used for NGC 4472, 
scaled to a smaller r^- We found that the amount 
of mass dropout M^g is higher in smaller ellipticals 
relative to the total baryonic mass M,. This is 
opposite to the trend observed in the fundamental 
plane. However, if hot gas and dark matter in the 
outer halos of smaller ellipticals is tidally stripped 
in group environments, as suggested by Mathews 
& Brighenti (1998b), then the mass of cooled gas 
would be reduced and its effect on the fundamen- 
tal plane would be reduced. Nevertheless, expla- 
nations of the deviations of the fundamental plane 
from virial scaling must recognize the possible ad- 
ditional influence of dropout mass, regardless of 
the trend of Mcg/M* with M*. 

6.6. Conclusions 

Using simple spherical gas-dynamical models 
for the evolution of interstellar gas and data from 
the well-observed elliptical NGC 4472, we reach 
the following conclusions: 

(1) If the hot interstellar gas cools only in the 
very center of NGC 4472 for ~ 10 Gyrs, the to- 
tal accumulated mass there would be > 10 times 
larger than the mass of the central black hole ob- 
served. If such large concentrated masses were 
generally present in bright ellipticals, interstellar 
gas in r ^ 0.1 re would be compressed and heated 
to temperatures > 1 keV. Such hot thermal cores 
are not generally observed. We conclude that the 
cooling dropout in massive ellipticals must occur 
before the gas reaches the galactic center. The 
hypotheses of distributed mass dropout and low 
mass star formation were proposed many years 
ago (Fabian, Nulsen & Canizares 1982; Thomas 
1986; Cowie & Binney 1988; Vedder, Trester, & 
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Canizarcs 1988; Sarazin & Ashe 1989; Fcrland, 
Fabian, & Johnstone 1994). However, these his- 
torical arguments were generally based on the no- 
tion that mass dropout would help reduce com- 
puted X-ray surface brightness profiles at small 
projected radii, as required by the observations, 
but wc have shown here that in some cases en- 
hanced dropout at small galactic radii can cause 
to increase, not decrease. (The total bolomet- 
ric X-ray luminosity should always be lower 
in distributed cooling models since the hot gas 
experiences only a fraction of the galactic poten- 
tial.) The best arguments for distributed cooling 
dropout are (i) limits on the central black hole 
mass and (ii) the absence of rotational flattening 
in X-ray images. 

(2) We have considered a wide variety of possible 
mass profiles for the radial deposition of cooled in- 
terstellar gas in NGC 4472. The dropout mass is 
assumed to be optically dark, consistent with the 
formation of very low mass stars. In every case the 
(stellar plus dropout) mass to light ratio at ^ r,,/3 
significantly exceeds the mass to light ratio deter- 
mined from stellar velocities. If the dropout mass 
is optically dark, dynamical mass to light ratios 
in luminous ellipticals should be substantially en- 
hanced by dark baryonic matter. In this case the 
true M/Lb for luminous stars may be ~ 30 per- 
cent smaller than published values. 

(3) The excellent agreement between the X-ray 
and "stellar" mass in NGC 4472 shown in Figure 1 
(and also NGC 4649) in the range O.lrg ^ r <^ Ir^ 
may be a coincidence if our estimates of the cool- 
ing flow dropout mass are correct and if this mass 
is non-luminous. 

(4) Dynamical mass to light determinations within 
Ve refer to a superposition of two stellar popula- 
tions: an old luminous population with a dc Vau- 
couleurs galactic profile and a younger population 
having a bottom-heavy IMF and a different galac- 
tic mass profile. If the younger population is opti- 
cally dark, the mass to light ratio is not likely to 
be constant with galactic radius within rg- 

(5) It is not possible to reduce the total amount 
of mass deposited from the cooling flow simply 
by lowering the specific rate of stellar mass loss 
a*(t) = {dMtf/dt)/M^. For all reasonable power 
law initial mass functions we show that a* oc 
[M^l Lb)~^ ■ For given Lb, larger stellar masses 
M*(r) must accompany lower values of a* so 



the total amount of mass ejected from stars ex 
/ at,{t)M^dt is nearly independent of the IMF. 

(6) Among the models we consider, those with 
centrally concentrated mass dropout perform best 
in minimizing the overall disagreement with the 
central M/Lb and the X-ray determined mass in 
O.lre r ^ Ire- Constant q models in which the 
dropout is proportional to the local gas emissivity 
at every radius deposit less mass in O.lre ^ r < 
Ire, but may have gas temperatures that are too 
low {q > 4) or central masses that are too large 
(9^1). 

(7) Even in the presence of mass dropout, the com- 
puted central interstellar gas densities and X-ray 
surface brightnesses {R) are generally too large. 
Such deviations would be reduced if the hot gas 
is partially supported in r < O.lre by magnetic 
or other non-thermal pressure associated with the 
extended radio source. 

(8) If the stellar population formed from cooled 
interstellar gas extends to intermediate masses, 
~ 1 — 2 Mq, its mass to light ratio may blend 
with that of the older population. In this case 
the dropout mass would already be represented in 
the de Vaucouleurs profile representing the stellar 
mass distribution in our models. Dynamical deter- 
minations of M/Lb would be a weighted mean of 
the two populations. If the dropout stellar popu- 
lation is luminous, some of the difficulties we have 
discussed here would be alleviated, but not those 
regarding T.x{R)- 

(9) If the infiuence of rotation and non-thermal 
pressure can be understood, high resolution im- 
ages of the central regions of elliptical galaxies us- 
ing the Chandra (AXAF) satellite may detect the 
presence of dark baryonic dropout material and an 
accurate determination of the mass to light ratio 
of the old stellar population. 
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Table 1 

MASS DROPOUT MODELS FOR NGC 4472*^ 



Model 


9o 




m 










"tot I \t 


A„/i(^)* 


Meent" 


r ros m 








(kpc) 




(10" Mo) 


(10" Mo) 


(10" Mq) 


(10" Mq) 


(Mq/Lbq) 




(10" Mo) 


(10*° erg s"^) 


(10*° erg s 


1 









4.64 


4.64 


1.35'' 


1.36 


13.70 


0.00 


4.64 


28.6 


28.6 


2 









4.66 


4.66 


1.81 


1.70 


13.71 


0.34 


4.66 


21.7 


21.7 


3 


1 






36.1' 


2.64 


1.61 


1.26 


12.19 


0.19 


0.254 


18.7 


36.4 


4 


4 






81. ij 


1.36 


1.48 


0.97 


11.21 


0.10 


0.013 


8.8 


42.6 


5 


4 


15 


1 


4.65 


1.55 


1.50 


1.01 


11.36 


0.11 


0.013 


13.9 


27.7 


6 


4 


2 


1 


4.65 


3.76 


1.72 


1.51 


13.03 


0.27 


0.015 


22.1 


33.5 


7 


4 


0.8 


2 


4.65 


4.65 


1.81 


1.76 


13.71 


0.34 


0.017 


29.1 


43.3 


8' 


4 


2 


1 


4.60 


3.65 


1.29 


1.17 


9.74 


0.40 


0.012 


20.0 


28.1 



^All masses evaluated at time tn = 13 Gyrs. 

^Total cooled mass within 1 Mpc. 

*^Mass of cooled gas within re/3 = 2.86 kpc. 

Total mass within re/3; M*,d^(r3/3) = 1.36 X 10^^ Mq. 
■^Mass of cooled gas at r^/S evaluated using hydrostatic equilibrium. 
^M*, dh/ Lb — 10.18 is the value for the stars and dark halo for Models 1-7. 
^Relative contribution of dropout mass to M/Lb at re/3. 

Mass cooled into central gridzone; the central black hole in NGC 4472 has mass Mhh = 0.29 X 10" Mo . 
'Mass within 100 Icpc is 4.72 X lo" Mq 
'Mass within 100 Icpc is 4.96 X lO"^" Mq 
''Mass of cooled gas is not included, 
'with lower M/Lb for old stars. 

"ROSAT X-ray luminosity (0.2 - 2 keV) of background flow 
^Total ROSAT X-ray luminosity including emission from dropout 



Table 1 

MASS DROPOUT MODELS FOR 

NGC 4472 - TABLE 1 
CONTINUED HERE BECAUSE OF 
AASTEX BUG 



Model 


r ros ni 


T ros n 
^x,tot 




(lO^o erg s-i) 


(lO^o erg s-i) 


1 


28.6 


28.6 


2 


21.7 


21.7 


3 


18.7 


36.4 


4 


8.8 


42.6 


5 


13.9 


27.7 


6 


22.1 


33.5 


7 


29.1 


43.3 


8 


20.0 


28.1 
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Fig. 1. — (a) Interstellar gas temperature and den- 
sity in NGC 4472. top: fit to ROSAT data for 
r(i?); bottom: fit to ROSAT data for n(r) (open 
eireles) and Einstein data, (closed eircles); (b) solid 
curve: Mx{r) determined from above data using 
equation (1); long-dashed curve: de Vaucouleurs 
mass profile M*(r) assuming M^^/Lb = 9.2. 

Fig. 2. — (a) Comparison of observed hot gas den- 
sity and temperature in NGC 4472 with various 
hydrodynamical models at time = 13 Gyrs. 
The source of observed data in all panels is iden- 
tical to that in Figure 1. Light solid lines: density 
and temperature of the background flow as func- 
tions of physical radius r; heavy solid lines: appar- 
ent density and temperature as functions of r in- 
cluding emission from cooling regions; dotted lines: 
apparent gas temperature as function of projected 
radius, similar to the temperature data points, (b) 
Comparison of ROSAT and Einstein X-ray surface 
brightness distributions in NGC 4472 with various 
models. Light solid lines: emission just from the 
background cooling flow; heavy solid lines: com- 
bined emission from background and distributed 
cooling regions. 



computed models using equation (1); dashed lines: 
actual mass profiles Mtot{r) for each model; dot- 
ted curves: de Vaucouleurs profiles M,(r) based 
on imiform mass to light ratios, M/Lb ~ 9.2 for 
Model 6 and M/Lb = 6 for Model 8. Circles: 
data points show the apparent mass distribution 
Ma,(r) using NGC 4472 observations and equation 
(!)• 

Fig. 6. — Plot of the specific mass loss rate a*(t„) 
(in sec~^) against the mass to light ratio M^JLb in 
solar units for 84 power law initial mass functions 
with varying slopes and mass cutoffs as described 
in the text. 



Fig. 3. — Baryonic and non-baryonic mass dis- 
tributions in NGC 4472. Light solid line: de 
Vaucouleurs profile of dynamical mass based on 
M/Lb = 9.2; heavy solid line: NFW dark halo 
profile for best fit to NGC 4472 data; short dashed 
line: cooled dropout mass Meg for Model 1; long 
dashed line: cooled dropout mass M^g for Model 4; 
dotted line: cooled dropout mass Meg for Model 6; 
short dash-dotted line: cooled dropout mass Meg 
for Model 7; long dash-dotted line: cooled dropout 
mass Meg for Model 8. 



Fig. 4. — (a) Comparison of observed hot gas den- 
sity and temperature in NGC 4472 with Model 8 
at time f„ = 13 Gyrs. Data and curve designations 
are identical to those in Figure 2a. (b) X-ray sur- 
face brightness profiles at i„ = 13 Gyrs for Model 
8. Data and curve designations are identical to 
those in Figure 2b. 



Fig. 5. — Mass profiles at tn = 13 Gyrs for Model 
6 (upper panel) and Model 8 (lower panel). Solid 
lines: X-ray mass profiles Mx{r) derived from the 
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